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We present detailed investigations of the magnetic properties of an Fe monolayer on W and Ta 
(110) surfaces based on the ab initio screened Korringa-Kohn-Rostoker method. By calculating 
tensorial exchange coupling coefficients, the ground states of the systems are determined using 
atomistic spin dynamics simulations. Different types of ground states are found in the systems as a 
function of relaxation of the Fe layer. In case of W(llO) substrate this is reflected in a reorientation of 
the easy axis from in-plane to out-of-plane. For Ta(llO) a switching appears from the ferromagnetic 
state to a cycloidal spin spiral state, then to another spin spiral state with a larger wave vector and, 
for large relaxations, a rotation of the normal vector of the spin spiral is found. Classical Monte 
Carlo simulations indicate temperature-induced transitions between the different magnetic phases 
observed in the Fe/Ta(110) system. These phase transitions are analyzed both quantitatively and 
qualitatively by finite-temperature spin wave theory. 
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I. INTRODUCTION 


The Dzyaloshinsky-Moriya interactioiP^^ 
between local magnetic moments has a great 
impact in spintronics applications through the 
formation of chiral spin structures like ma^ 
netic skyrmion^^ and chiral domain wallspl 
while it may also lead to an asymmetry in 
the magnon spectrum of ferromagnetic thin 
films, as was shown theoreticall}^^ and ex¬ 
amined in spin-polarized electron energy loss 
spectroscopy experiments^ for Fe/W(110). 
Spin-polarized scanning tunneling microscopy 
experiments enabled the real-space observa¬ 
tion of spin spiral orderings at low temper¬ 
atures in several ultrathin films such as Mn 
monolaw on W(110),^Pd/Fe double-layer on 
Ir(lll),fil Cr monolayer on W(110)P^ and Fe 
double-layer on W(llO).^^^^ 


A double-layer of Fe on W(llO) shows un¬ 
usual phase transitions when the temperature 
is increased. While the monolayer is ferromag¬ 
netic up to Tc ~ 230 in the double-layer 
the low-temperature spin spiral phase disap¬ 
pears at around 200 developing an in¬ 
plane ferromagnetic state as in the case of the 
monolayer, which persists up to ~ 450 K.^ 
This is in agreement with the asymmetry of 
the spin wave spectrum found in Ref. [7] 
at T ^ 300 K, since the spectrum around a 
cycloidal spin spiral ground state would be 
symmetric if the Dzyaloshinsky-Moriya in¬ 
teraction were perpendicular to the plane of 
the spiral.ll^ll^ Using the experimentally ob¬ 
tained wavelength of the low-temperature spi¬ 
ral state it was possible to find micromag- 
netic exchange (spin stiffness), Dzyaloshinsky- 
Moriya and anisotropy para meters describ¬ 
ing this type of order.^^^^^ However, both 



2 


micromagneticP^^ and atomisticP^ ab initio 
calculations indicated a ferromagnetic ground 
state in the system. For an Fe monolayer 
on W(llO), theoretical calculation^^^^f^ agree 
with experiments^ in determining an in¬ 
plane ferromagnetic ground state. For a Mn 
monolayer on W(llO), Ref. [9] provided con¬ 
sistent experimental and theoretical descrip¬ 
tions of the spiral ground state. Ab initio 
calculations^ and experimentS^ are also in 
agreement about the spiral ground state of Cr 
monolayer on W(llO). 

Various types of magnetic ground state con¬ 
figurations were found by ab initio calculations 
in an Fe monolayer on the (100) surface of 
Wi-a^Ta^^ (0 < X < 1) alloyS^^^ as a func¬ 
tion of Ta concentration x, ranging from an 
antiferromagnetic state on pure W to a ferro¬ 
magnetic state on pure Ta. Both W and Ta 
have bcc lattice structure but the lattice con¬ 
stant of Ta is about 4.3% larger than that of W 
(<^Ta = 3.301 A and aw = 3.165 A ). This dif¬ 
ference was taken into account by calculating 
the lattice constant of the alloy, but the relax¬ 
ation of the Fe layer towards the top substrate 
layer was kept fixed during the calculations at 
the value determined for Fe/W(100), although 
Fe should have a larger inward relaxation in 
the case of Ta with the larger lattice constant. 
For different relaxations, Fe on Ta(lOO) may 
have either ferromagnetic or antiferromagnetic 
ground state as shown in Ref. m- 

In this paper we examine the magnetic 
ground state of an Fe monolayer on W and Ta 
(110) surfaces as a function of the relaxation 
of the Fe layer with respect to the top sub¬ 
strate layer. The electronic structure calcula¬ 
tions were performed by using the relativistic 
screened Korringa-Kohn-Rostoker method.!^ 
For the determination of the magnetic ground 
state we mapped the spin system onto a 
generalized classical Heisenberg model, where 
the parameters are taken from the relativis¬ 
tic generalizatiorP^ of the method of infinites¬ 
imal rotations introduced by Liechtenstein et 
al\M The ground state of the system was found 
by atomistic spin dynamics si mulat ions based 
on the Landau-Lifshitz-GilberfP^^^ equations. 
These results are described in Sec. im 

Besides changing the relaxation, thermal 
fluctuations may also induce transitions be¬ 


tween the different types of ordered states 
found in these systems. Classical Monte Carlo 
simulations were performed using the previ¬ 
ously obtained spin model to find these tran¬ 
sitions. In one of the transitions found in Fe 
monolayer on Ta(llO) the increasing temper¬ 
ature drives the system from the ferromag¬ 
netic ground state into a non-collinear spin 
spiral state. Most likely, this transition is 
driven by the Dzyaloshinsky-Moriya interac¬ 
tions and the easy-axis anisotropy in the sys¬ 
tem. Such a transition was already stud¬ 
ied in Refs. [33] and [34] using a Ginzburg- 
Landau model, which is, however, unsuitable 
for employing Heisenberg model parameters 
obtained from ab initio calculations. 


Instead of relying on a continuum model, 
we used spin wave expansion to describe 
the transition between the different ordered 
states. Th is me thod was found to be a 
powerful too R^ I ^^ I for explaining a transition 
from a low-temperature ferromagnetic to a 
high-temperature helical state in bulk Dy. 
In the present work we incorporated the 
Dzyaloshinsky-Moriya interaction into such 
an analysis, which was unnecessary in bulk 
systems with an inversion center, but it 
plays an important role in case of ultrathin 
films. We also used the spin wave expan¬ 
sion technique to handle higher order terms 
(magnon-magnon interactions) perturbatively, 
since perturbation theory makes it possible to 
estimate the temperature where the system 
reaches the paramagnetic state. This method 
was originally used to calculate the Curie tem¬ 
perature in a simple cubic lattice described by 
a ferromagnetic Heisenberg model.l^ By us¬ 
ing a simplified model Hamiltonian consistent 
with the different types of ground states fou nd 
in an Fe monolayer on Ta(llO), in Sec. HI we 


present a detailed analysis of the temperature- 
induced magnetic phase transitions and relate 
the results to those obtained from Monte Carlo 
simulations. 
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II. MAGNETIC STATES AND PHASE 
TRANSITIONS IN AN FE 
MONOLAYER ON W AND TA (110) 
SURFACES 

A. Ab initio calculation of collinear 
magnetic states 

For the ab initio calculations we used the 
relativi stic scre ened Korringa-Kohn-Rostoker 
method p^ * ^^ * ^^ * using the local spin density ap¬ 
proximation and the atomic sphere approxi¬ 
mation. First we performed calculations for 
W and Ta bulk with the lattice constants 
aw = 3.165 A and axa = 3.301 A, respec¬ 
tively. The layered systems considered for the 
deposited Fe monolayers comprised eight lay¬ 
ers of bulk atoms, one layer of Fe and three 
layers of empty spheres, sandwiched between 
the semi-infinite bulk calculated in the previ¬ 
ous step and a semi-infinite vacuum. Theo¬ 
retical calculations using the full-potential lin¬ 
earized augmented plane-wave method give re¬ 
laxation values betwee n 12 — 13% for an Fe 
monolayer on W(110),I^^®®1 while the ex¬ 
perimental values are in the range of 7 — 
13^[43}{45] Ta(llO) Fe should have an even 
larger relaxation due to the larger lattice con¬ 
stant. Therefore the calculations were per¬ 
formed for different values of the distance be¬ 
tween the Fe monolayer and the top bulk 
monolayer, adjusting the Wigner-Seitz radius 
of the atomic spheres related to the Fe atoms 
correspondingly. Both for W and Ta, the 
relative relaxation with respect to the ideal 
distance between bee(110) atomic layers was 
changed between 10% and 17%. All the 
atomic layers but the Fe layer were kept at the 
ideal lattice geometry since calculationJ^^®!] 
indicate that the W-W relaxations are be¬ 
low 1% even between the topmost W mono- 
layers. We determined the potential and 
the exchange-correlation magnetic field self- 
consistently, serving as an input to the evalua¬ 
tion of the exchange coefficients, see Sec. El 
The spin and orbital magnetic moments 
obtained from the ab initio calculations are 
listed in Table The sum of the spin and or¬ 
bital moments in the Fe layer on W(llO) for 
13% inward relaxation compares within 10% 
to the total magnetic moments given in the 


literature.The induced moments in the 
topmost W layer are antiparallel to the Fe mo¬ 
ments, in agreement with Refs. [40] and |41| . 
but they are parallel in the next two W lay¬ 
ers. It is worth noting that the spin and or¬ 
bital moments are parallel for the W atoms 
although the W d-shell is less than half-filled, 
which indicates a violation of Hund’s third 
rule, cf. Ref. [41]. Apparently, this is not 
the case for Ta. It is also notable that the 
induced moments of the Ta atoms are larger 
than those of the corresponding W atoms. 
Ref. [42] agrees with our calculation inasmuch 
as increasing the relaxation decreases the mag¬ 
netic moments of the Fe atoms, most likely due 
to the increased hybridisation between the Fe 
and the substrate layers. 

B. Calculated exchange interactions 

Using the self-consistent potentials obtained 
before, the relativistic torque method^ was 
employed to map the energy of the magnetic 
system onto a generalized Heisenberg model, 

^ ^ E , (1) 

i,j i 

where i, j and a, (3 label lattice sites and Carte¬ 
sian indices, respectively, Sf are the compo¬ 
nents of the unit vector representing the ori¬ 
entation of the spin at lattice site i, while 

and stand for the matrix elements of the 
exchange coupling tensors and of the second- 
order on-site anisotropy energy tensors. The 
relativistic torque method relies on the mag¬ 
netic force theorem and requires the calcu¬ 
lation of coupling coefficients around differ¬ 
ent collinear reference states for at least three 
linearly independent magnetization directions, 
since for a given direction, only those compo¬ 
nents of the Jij tensors can be obtained which 
lie in the plane perpendicular to the magneti¬ 
zation. In particular, we considered the mag¬ 
netization directions [110], [001], and [110]. 
The spins in a given layer must be ferromag- 
netically aligned, but the antiferromagnetic or¬ 
dering between the different layers was taken 
into account. To perform the necessary inte¬ 
grations, 16 energy points were taken along a 
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Fe/W(110) 

spin moment (/xb) orbital moment (/xb) 


relaxation 

Fel 

W1 

W2 

W3 

Fel 

W1 

W2 

W3 

10% 

2.355 

-0.164 

0.007 

0.003 

0.180 

-0.027 

0.001 

-0.001 

13% 

2.244 

-0.164 

0.012 

0.003 

0.169 

-0.018 

0.005 

0.000 

15% 

2.181 

-0.161 

0.017 

0.004 

0.162 

-0.014 

0.007 

0.001 

17% 

2.122 

-0.156 

0.022 

0.004 

0.156 

-0.011 

0.010 

0.002 


Fe/Ta(110) 

spin moment (/xb) orbital moment (/xb) 


relaxation 

Fel 

Tal 

Ta2 

Ta3 

Fel 

Tal 

Ta2 

Ta3 

10% 

2.587 

-0.278 

-0.037 

-0.027 

0.100 

0.031 

0.005 

0.003 

13% 

2.520 

-0.310 

-0.045 

-0.030 

0.097 

0.036 

0.006 

0.003 

15% 

2.466 

-0.333 

-0.046 

-0.027 

0.094 

0.040 

0.006 

0.002 

17% 

2.406 

-0.358 

-0.044 

-0.023 

0.090 

0.044 

0.006 

0.001 


TABLE L Calculated spin and orbital moments in the Fe layer and in the topmost three substrate layers 
of W(llO) and Ta(llO) surfaces for selected values of relaxations of the Fe layer. 


semicircle contour in the upper complex semi¬ 
plane, and from 204 up to 6653 /c-points were 
sampled in the Brillouin zone, gradually in¬ 
creasing for energies approaching the Fermi 
level. 

The isotropic part of the exchange tensors 
between the Fe atoms, 

= ( 2 ) 

a 

is shown in Fig. for W and Ta surfaces and 
different relaxation values. Note that with the 
sign convention of Eq. Q, Jij < 0 and Jij > 0 
indicate ferromagnetic and antiferromagnetic 
couplings, respectively. In case of W(llO) sur¬ 
face, the nearest-neighbor ferromagnetic cou¬ 
pling is fairly insensitive to the relaxation, 
while the next-nearest-neighbor coupling (at 
the distance of one lattice constant) is anti¬ 
ferromagnetic for lower relaxations, but be¬ 
comes ferromagnetic above 15% relaxation. 
For Ta(llO) surface, the exchange couplings 
for the two nearest neighbors are ferromag¬ 
netic for all considered values of relaxations. 
The weaker nearest-neighbor interaction de¬ 
creases and the next-nearest-neighbor inter¬ 
action increases in size with increasing relax¬ 
ation. Also notable is the increasingly anti¬ 
ferromagnetic character of some further (third 
and fifth) neighbor couplings with increas¬ 
ing relaxation, which will give rise to the 
formation of a short wavelength spin spiral 




distance (a-^^) 


FIG. 1. (color online) Calculated isotropic ex¬ 
change interactions Jij obtained from the rela¬ 
tivistic torque method, for (a) W(llO) and (b) 
Ta(llO) surfaces and different values of relaxations 
of the Fe layer. 
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FIG. 2. (color online) Calculated components of the Dzyaloshinsky-Moriya vectors, Dfj and , in an 
Fe monolayer on (a)-(b) W(llO) and (c)-(d) Ta(llO) surfaces for different values of relaxations of the Fe 
layer. Displayed are the values for the pairs with Rfj > 0 and > 0, while the Dzyaloshinsky-Moriya 
vectors for the rest of the pairs can be obtained from the symmetry relations described in the text. 


along the [110] direction in Fe/Ta(110), see 
Sec. |I^ In particular, this might happen 
since the strong ferromagnetic coupling be¬ 
tween the next-nearest neighbors does not play 
a role in the formation of the spiral state since 
it only couples spins along the [ 001 ] direction 
(see coupling J 2 in Fig. [^. 

The antisymmetric part of the exchange 
tensors between the Fe atoms is shown in 
Fig. |2] in terms of the components of the 
Dzyaloshinsky-Moriya vectors, 

= (3) 

/ 3,7 

AccordiM to the symmetry rules set up by 
Moriya,^ all the Dzyaloshinsky-Moriya vec¬ 
tors lie in the (110) plane. Note that the x 
and y directions correspond to the [ 110 ] (long) 
axis and to the [ 001 ] (short) axis, respec¬ 


tively. The components of the Dzyaloshinsky- 
Moriya vectors are only drawn for neighbors 
with Rfj > 0 and R^^ > 0. The compo¬ 
nents for the related neighbors can be obtained 
by symmetry: for 

{Dfj,-Dy^) for i-Rfj,Rl) and (-Df.,-Z)f.) 
for {—Rfj, —Rij)- Dfj is, therefore, only finite 
between atoms which have a finite distance 
along the [ 001 ] (y) direction; for example, the 
atoms at ^/2a distance are located along the 
[110] {x) axis, thus Dfj = 0. Similarly, 
is only finite if Rfj 7 ^ 0. The Dzyaloshinsky- 
Moriya interactions are comparable in magni¬ 
tude to the isotropic exchange interactions and 
they also show oscillating behavior. 

The presence of the Dzyaloshinsky-Moriya 
interactions may stabilize spin spiral states 
and the sign of the components of the 
Dzyaloshinsky-Moriya vectors determines the 
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chirality of the spin spiral. Let q be the wave 
vector of the spiral, n the normal vector of the 
monolayer pointing outwards from the sub¬ 
strate, and introduce the vector % = x Sj 
such that q{Rj — Ri) > 0 , where Ri and Rj 
are the position vectors of neighboring spins 
in the lattice. Note that for cycloidal spin spi¬ 
rals the direction of x is independent of the 
choice of the lattice sites i and j. Follow¬ 
ing Refs. m and m, a cycloidal spin spi¬ 
ral is called right-rotating when the vectors 
(gr, X, n) form a right-handed system. If they 
form a left-handed system, the spin spiral is 
called left-rotating. With our sign conven¬ 
tion and only taking into account the largest 
Dzyaloshinsky-Moriya interactions in both di¬ 
rections, in the case of W substrate the Dfj 
component prefers a right-rotating spiral along 
the [ 001 ] direction and the D^- component 
prefers a left-rotating spiral along the [110] di¬ 
rection. This is in agreement with the results 
in Ref. [19] and the chirality of the spin spiral 
state along the [001] direction in double-layer 
Fe on W(llO).^ For Ta substrate, the sign 
of the largest D^- vector component is flipped 
compared to the case of W substrate. This 
means that the Dzyaloshinsky-Moriya interac¬ 
tions prefer left-rotating spirals in an Fe mono- 
layer on Ta(llO) along both the [001] and [110] 
directions. 


C. Ground states obtained from spin 
dynamics simulations 


After obtaining the coupling coefficients 
from collinear configurations, we performed 
atomistic spin dynamics simulations to find 
the ground states of the systems. These are 
based on the numerical solution of the Landau- 
Lifshitz-Gilbert equations. 


dtSi = --f'Mi-a^'SiXMi, (4) 

with 7 ' = T^^ 2 m gyromagnetic coeffi¬ 
cient {g the g-factor, e the magnitude of charge 
and m the mass of the electron) and a the 
dimensionless Gilbert damping factor. The 
torque Mi acting on the spin vector Si is de¬ 
fined as 


Mi = SiX 


1 dH 

rrii dSi 


( 5 ) 


and rrii is the magnitude of the magnetic mo¬ 
ment of the atom at site i, associated with 
the spin magnetic moment from the ab initio 
calculations in Sec. EE while H is the spin 
Hamiltonian in Eq. W 

We also calculated the exchange couplings 
between the Fe atoms and the atoms in the 
topmost bulk layer which had the largest in¬ 
duced moment, see TableHowever, we found 
that including these couplings did not change 
the ground state considerably, they just give 
rise to an antiparallel alignment of the induced 
moments with respect to the neighboring Fe 
moments. This implies that for the considered 
systems only the stable Fe moments are rele¬ 
vant to be included into the Hamiltonian 0 . 
This feature is essential since the quasiclassi- 
cal description 0-0 is shown to be a reliable 
description for the rigid moments,!^ but it is 
probably not valid for the induced moments. 

Starting the spin dynamics simulations from 
a random initial configuration, the system will 
generally converge to a metastable equilibrium 
state, that is to a local energy minimum. How¬ 
ever, this configuration may not be the ground 
state - the global energy minimum -, there¬ 
fore the determination of the ground state may 
require multiple runs. It was found that a 
random initial state often leads to a spin spi¬ 
ral state, even if it has slightly higher energy 
than the ferromagnetic state. Furthermore, 
the obtained equilibrium states may contain 
skyrmion-like local excitations which are sta¬ 
ble with respect to the dynamics of the system, 
but represent a positive energy correction com¬ 
pared to the ground state. 

The energies obtained from the spin dynam¬ 
ics simulations with the Hamiltonian 0 are 
shown in Fig. |^a) as a function of the relax¬ 
ation of the Fe layer in case of W(llO). The 
ground state energy of the system is compared 
to the energies of the ferromagnetic align¬ 
ments along the main crystallographic direc¬ 
tions [110], [001] and [110]. The ground state 
of the Fe monolayer on W(llO) was found to 
be ferromagnetic for all relaxations, however, a 
reorientation transition occurs at around 15% 
relaxation of the Fe layer from the in-plane 
[ 110 ] direction to the out-of-plane [ 110 ] direc¬ 
tion. The in-plane easy axis at the experimen¬ 
tally observed relaxation value 13% is in agree- 
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10% 11% 12% 13% 14% 15% 16% 17% 


relaxation 



relaxation 


FIG. 3. (color online) Energies per Fe spin of an Fe 
monolayer (a) on W(llO) and (b) on Ta(llO) for 
different magnetic states as a function of the relax¬ 
ation of the Fe layer obtained from spin dynamics 
simulations for a system consisting of iV = 64 x 64 
atoms with periodic boundary conditions. The en¬ 
ergy of the ground state (GS) is highlighted by 
blue solid line and the types of the ground state 
magnetic orderings are displayed for the whole 
range of relaxations. For the explanation of the 
different spin spiral states (SS I, SS II, SS III) see 
the text. 

merit with the experiments.^^ It is worth not¬ 
ing that a double-layer of Fe on W(llO) has 
an out-of-plane easy axispEH similarly to the 
case here for large relaxation. 

In Fig. I^b) the energies of the ferromag¬ 
netic states and also of different spin spiral 
states are shown for the Fe monolayer on 
Ta(llO). The energies of the cycloidal spiral 
states SS I and SS II were calculated in the 
homogeneous left-rotating spin spiral configu¬ 
ration, 

Si = (- sin {qRi ), 0 , cos {qRi )), ( 6 ) 


where the different spin components cor¬ 
respond to the directions (x, z) = 
([110], [001], [110]). The normal vector 
and rotational sense of the spirals chosen 
in Eq. ^ are consistent with the obtained 
ground states shown in Fig. The spiral 
energies were calculated for q values in the 
whole Brillouin zone, but only the q vectors 
along the [ 110 ] direction, denoted by 
showed complex behavior, see Fig. In 
Fig.H the energy difference between the spin 
spiral states and the ferromagnetic state along 
the [ 110 ] direction does not go to 0 as g ^ 0 
due to the anisotropy in the system. 

Fig. I^b) indicates phase transitions at 
10.5% relaxation from the ferromagnetic state 
with out-of-plane easy axis (FM) to the SS I 
spin spiral state, at 13.8% relaxation between 
the SS I and SS II states, and at 14.5% re¬ 
laxation between the SS II and SS III states. 
All the spiral states have a wave vector par¬ 
allel to the [ 110 ] direction, and all the spins 
in the spiral are confined to a plane. For 
the SS I and the SS II states, the spins are 
located in the [ 110 ] — [ 110 ] plane, forming a 
left-rotating cycloidal ^in spiral as in a Mn 
monolayer on W(llO),^ although it is clear 
from Fig. [^b) that the [110] direction is the 
hard axis since the ferromagnetic state along 
this direction has the highest energy. The 
plane of the spiral is thus clearly a consequence 
of the Dzyaloshinsky-Moriya interaction in 
the system which prefers spin spiral states 
oriented perpendicular to the Dzyaloshinsky- 
Moriya vector. For a spin spiral along the 
[ 110 ] direction, only the [ 001 ] component of 
the Dzyaloshinsky-Moriya interaction plays a 
role in the ground state energy, leading to the 
cycloidal spiral state resembling a Neel domain 
wall. 

The SS I state has a small wave number, the 
value of which increases continuously with in¬ 
creasing relaxation (see the squares in Fig. 
but jumps to the much larger wave number 
of the SS II spin spiral at relaxation 13.8%. 
The presence of spin spiral energy minima 
at different wave vectors and the transition 
between these minima is a consequence of 
the frustrated isotropic exchange interactions 
around these relaxations, see Fig. Bb). The 
SS III state has similar wave number to the 
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(a) 

FM 

f[ooi] 


[no] [no] 


[ 001 ] [ 110 ] 






1110 ] 


[ 001 ] [ 110 ] 

(d) 

SS III 

^A|001] 


[ 110 ] [ 110 ] 




[ 110 ] 


[001] [no] 


FIG. 4. (color online) Ground state magnetic configurations of an Fe monolayer on Ta(llO) for different 
relaxations: (a) FM state for 10.3%, (b) SS I state for 13.5% (A 5.98 nm), (c) SS II state for 13.8% 
(A ~ 0.83nm) and (d) SS III state for 15% (A ~ 0.81 nm). The SS I and SS II states differ in the 
wavelength of the spin spiral, while the SS II and SS III states mainly differ in the normal vector of the 
spiral. 



FIG. 5. (color online) Spin spiral energies per spin 
relative to the energy of the ferromagnetic state, 
calculated from the Heisenberg model parameters 
in the spin spiral configuration Eq. ^ for wave 
vectors along the [110] direction, Qx (given in units 
of ). The points at which the spin spiral ener¬ 

gies are calculated in Fig.are denoted by squares 
for SS I and circles for SS II. The inset shows a 
magnified view of the range 0 < qx < 0.15. 


SS II state, however, the anisotropy is strong 
enough to rotate the plane of the spiral out 
from the [110] — [110] plane, that is the normal 


vector [001] changes to a general direction in 
the [110] — [001] plane. The ground state ener¬ 
gies obtained from the spin dynamics simula¬ 
tions in Fig.j^b) are somewhat lower than the 
spin spiral energies presented in Fig. since 
due to the anisotropy the spiral can gain en¬ 
ergy by being deformed with respect to the 
perfect sinusoidal shape.l^^ This difference is 
the largest for the SS III state, but in that 
case this is also a consequence of the rotation 
of the normal vector of the spin spiral. 


D. Phase transitions at finite 
temperature using Monte Carlo 
simulations 

We examined the phase transitions in the 
systems also for fixed relaxations as a function 
of temperature, using classical Monte Carlo 
simulations with Metropolis dynamics. These 
phase transitions were expected to occur for 
relaxation values close to the transition points. 
The order parameter of the simulations was 
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FIG. 6. (color online) Temperature dependence of the order parameters, Eqs. of an Fe monolayer 

obtained from Monte Carlo simulations for lattice size A/" = 64 x 64. The wave number q is given in 
units of and is in all cases parallel the [110] axis, (a) W(llO) surface, 15% relaxation, ferromagnetic 
order parameter; (b) Ta(llO) surface, 10% relaxation; (c) Ta(llO) surface, 13.8% relaxation; (d) Ta(llO) 
surface, 15% relaxation, for q = 0.593750. 


defined 


with 


m. 


m 




m. 


(9), 


a=x,y,z 


:(g) 




( 7 ) 


(8) 


where ( ) denotes thermal average. As dis¬ 
cussed in Sec. II C[ the shape of the spiral 
state will differ from a perfect sinusoidal shape 
due to the anisotropy in the system. There¬ 
fore the order parameter for wave vector q 
does not perfectly fit this anharmonic spiral 
with the same wave vector due to the appear¬ 
ance of higher Fourier harmonics, but still it 
gives a good approximation to characterize the 

ordering .SSI 


The temperature dependence of the order 
parameters is shown in Fig. For the Fe 
monolayer on W(llO), see Fig. [^a), no re¬ 
orientation transition occurred in the system, 
although the relaxation value of 15% was 
close to the transition point. Similarly, no 
temperature-induced reorientation was found 
on the other side of the phase boundary, at 
15.2% relaxation. The paramagnetic state was 
reached at Tc ~ 350 K, somewhat higher than 
the experimentally determined critical tem¬ 
perature, Tc ~ 230 K.SSl 

In case of the Ta substrate several types 
of temperature-induced transitions happened 
between the different ordered phases before 
reaching the paramagnetic phase, if the cho¬ 
sen relaxation value was close to the phase 
boundaries shown in Fig. [^b). The SS I 
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phase turned out to be the most stable one 
against thermal fluctuations: systems with fer¬ 
romagnetic ground state at 10 % relaxation or 
with a SS II ground state at 13.8% relaxation 
turned into the SS I state, in both cases at 
around 130 K, as indicated by a change in 
the wave number of the order parameter in 
Fig. gb) and Fig. [^c), respectively. More¬ 
over, in case of the FM-SS I phase transi¬ 
tion a continuous increase of the wave num¬ 
ber can be inferred from Fig. [^b) above the 
critical temperature of the phase transition. 
For the case of a SS III ground state at 15% 
relaxation. Fig. |^d) shows that the com¬ 
ponent decreases with the temperature simi¬ 
larly to the order parameter m? in Figs. |^b)- 
(c). However, initially increases with the 
temperature, which is accompanied by a more 
pronounced decrease of m^. This indicates 
that the normal vector of the spin spiral ro¬ 
tates towards the y = [ 001 ] axis and at about 
80 K a phase transition to the SS II state oc¬ 
curs. The paramagnetic state was reached at 
Tc ^ 140 — 220 K in the case of Ta substrate 
depending on the relaxation. 


III. DESCRIPTION OF THE PHASE 
TRANSITIONS IN FE/TA(110) BASED 
ON SPIN WAVE EXPANSION 


In this Section, the temperature-induced 
phase transitions in the Fe monolayer on 
Ta(llO) surface will be discussed in terms of 
spin wave expansion. Keeping the same global 
coordinate system as in Sec. IIC, {x^y^z) = 
([ 110 ], [ 001 ], [ 110 ]), we will use a simplified 
model Hamiltonian, 

^ E + b E 


i,j i,j 


( 9 ) 



FIG. 7. (color online) Sketch of the lattice and 
the model parameters considered in Eq. for 
an Fe monolayer on Ta(llO). Jj denote effec¬ 
tive exchange couplings between the spin at site 
0 and its neighbors (see text). Equivalent neigh¬ 
bors are formed by mirroring on the xz and yz 
planes: there are four neighbors of types 1 and 
7, as well as two neighbors of types 2, 3 and 11. 
Only the nearest-neighbor Dzyaloshinsky-Moriya 
vector Di is taken into account, and it transforms 
as an axial vector. 


Since the spin spirals have a wave vector par¬ 
allel to the X axis, only such parameters are 
relevant which influence the ordering along 
this direction. These are the effective ex¬ 
change couplings denoted by Ji, J 2 ^ Jsj Jr, Jii 
and a Dzyaloshinsky-Moriya vector between 
the nearest neighbors Di parallel to the y 
axis (see Fig. [^. The isotropic couplings are 
summed up along the y axis: for example, J 3 
represents the coupling between the spin at 
site 0 and all the atoms which have the same 
X coordinate as the third neighbors. This is 
because the contributions of these Fe-Fe pairs 
add up in the energy of the spin spirals with 
wave vectors along the x axis. The anisotropy 
constants are chosen in agreement with the en¬ 
ergies of the ferromagnetic states along the dif¬ 
ferent axes in Fig. [^b). 

Within the spin wave expansion, the energy 
of the spin system is expanded around a stable 
equilibrium state using small spin deviations 
with respect to this state. To lowest order, 
the Hamiltonian can be written as 

Ho =Eo + J2 

k 


where Jij = Dij = ( 0 ,D^j, 0 ) with Dij > 


-D K 

J %•) Z 


'ij 

< 0 


0 for Rfj > 0 and Dij 
and Kx > 0 , that is 2 : is the easy axis and 
X is the hard axis. We choose the parame¬ 
ters such that the above Hamiltonian repro¬ 
duces the different phases found in Sec. |H C 


where Eq is the energy of the equilibrium state, 
the ak variables are the classical equivalents of 
bosonic spin wave annihilation operators and 
the spin wave energies, > 0 , stand for the 
energy corrections due to the spin excitations 
represented by a/.. 
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For Kx = Kz = 0^ a. homogeneous cy¬ 
cloidal spiral state in the xz plane with 
wave vector along the x axis, Si = 
(— sin (qQRi ), 0, cos (go-^i))? is either a stable 
or an unstable equilibrium state of the system. 
The energy per atom of the spin spiral is given 
by 

=ljiQo) - (11) 

with 

J(q)=^ (12) 

Dig) =J2 (13) 

The spin wave spectrum around a homoge¬ 
neous cycloid al sp iral state with wave vector 
qQ is given b> b^F^ I 

=\/C+{<l'^Qo)C- (q;qo), (14) 

with 

c+ (q; go) =\ [J (Q -<1o) + J(.Q + 9o)] 

[i-D (q + qo) - i£> (q - qo)] 

- J (qo) + i£> (qo), (15) 

C- (q; qo) =j (q) - j (qo) + i-D (qo), (16) 

where the excitations are indexed with the 
Fourier transformation wave vectors q. The 
equilibrium state is stable if both (7+ (q; Qo) 
and C- (q; Qq) non-negative for every q, 
which leads to real and non-negative spin wave 
frequencies.!^ The condition (7+ (q; Qo) ^ 
0 generally holds true if the wave vector 
Qo is close, but not necessarily equal, to 
the value for which Eq. 0 is minimized. 
Without Dzyaloshinsky-Moriya interactions, 
C- (q; Qo) >0 only holds if J (qo) is the 
global minimum of J (q). However, the pres¬ 
ence of the Dzyaloshinsky-Moriya interaction 
stabilizes several spiral states with different 
Qo values by achieving (7_ (q;Qo) ^ 0, even 
ones which do not minimize Eq. This 

leads to the appearance of metastable states 
for which the spin wave expansion ( p!Q| ) applies. 
The presence of the anisotropy may also sta¬ 
bilize these spiral states, either by introducing 


a hard axis perpendicular to the spiral plane 
{Ky > 0 in our model, cf. Ref. [35]) or by 
introducing an easy axis in the plane of the 
spiral {Kz < 0, cf. Ref. [36]b 

At finite temperatures, the free energy per 
atom of a system described by the spin-wave 
Hamiltonian (10) can be expressed as 


N ~ N ^ N 


■C(T), (17) 


where C (T) does not depend on the parame¬ 
ters of the equilibrium state Eq and This 
expression can describe a transition between 
two different stable equilibrium states spec¬ 
ified by parameters Eo^uJk and If 

Eq < Eq and the relation Incj/. > 
applies, then the system will switch from the 
first state to the second one at the temperature 


^B^trans — 


E'q-Eq 


Ek InWfc-EfclnWfe' 


(18) 


The quantum version of this method was ap¬ 
plied in Refs. [35] and [36] to describe the 
transition from a ferromagnetic to a spin spi¬ 
ral state in Dy. It should be noted that this 
method only gives numerically good transition 
temperatures if the temperature itself is small, 
since the spin wave expansion for the free en¬ 
ergy (17) becomes less accurate as the temper¬ 


ature is increased. 

A way of including a perturbative correction 
in the calculations is by writing the free energy 
as 


E —Eq + CJ/cU/c + - Pi 




k,k' 


—k^T ^ Inn/c, 


(19) 


where Pkk' is a symmetric matrix represent¬ 
ing higher order corrections to the energy ( 10 ) 


and Uk is the occupation number of the spin 
wave with energy uj^. Minimizing (19) with 
respect to rik leads to self-consistent equations 
which have real nonnegative solutions only for 
T < Tc, giving an estimate of the transi¬ 
tion temperature into the paramagnetic phase. 
This method was originally applied in Ref. m 
to find the Curie temperature of a Heisenberg 
ferromagnet on a simple cubic lattice. 
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A. The FM-SS I transition 


Based on the ab initio calculations, we 
chose different sets of model parameters which 
are close to the transition points, and em¬ 
ployed the spin wave expansion described 
above to obtain the possible phase transitions 
as a function of temperature. The calcula¬ 
tions were compared to Monte Carlo simu¬ 
lations using the Metropolis algorithm. For 
10 — 11% relaxations the spin spiral energy in 
Fig.|5] had a single minimum, which we re¬ 
produced by choosing a nearest-neighbor fer¬ 
romagnetic coupling Ji = — 2.0mRyd and 
a Dzyaloshinsky-Moriya parameter of Di = 
O.dmRyd. Kz = —0.22mRyd was used to 
move the energy of the ferromagnetic state 
below the minimum of the spin spiral energy 
curve. 

As pointed out in Sec. |II C| the equilibrium 
spin spiral states of the system are no longer 
perfect sinusoidal waves, since the spins will 
prefer the 2 : direction over the x axis. After 
finding this exact equilibrium state numeri¬ 
cally, the spin wave expansion was first per¬ 
formed on the basis of Eq. ( pT] ), that is for 
a system of free spin waves. The results are 
given in Table for a lattice size of N = 
128 X 64. The energies of the equilibrium states 
increase with increasing wave vector, with the 
ferromagnetic state {q = 0) being the ground 
state. Although the size of the lattice infiu- 
ences the allowed wave vectors in the Brillouin 
zone, the ground state rema ins fe rromagnetic 
even in the continuum limilP^^ correspond¬ 
ing to an infinite lattice, if the anisotropy is 
large enough. The free energy correction per 
spin due to free spin waves de¬ 

creases when the wave number of the spin spi¬ 
ral increases, leading to the expected transi¬ 
tion from the ferromagnetic to the spiral state 
with increasing temperature. After this transi¬ 
tion, the wave number of the equilibrium spin 
spiral gradually increases. This change is con¬ 
tinuous in the continuum limit, therefore the 
spiral orderings with different wave vectors do 
not actually represent different phases. 

Including perturbation corrections in the 
calculations on the basis of Eq. (19) makes 
it possible to give an approximation for Tc, 
where the equilibrium state loses its stabil- 


n f ^ 

A (nm) 

^ (mRyd) 

-^sw 

TSUs (K) 


Nk^T 

0.000000 

oc 

-4.2200 

1.9630 

0.0 

0.015625 

29.88 

-4.2176 

1.9531 

38.2 

0.031250 

14.94 

-4.2152 

1.9438 

40.7 

0.046875 

9.96 

-4.2124 

1.9353 

51.9 

0.062500 

7.47 

-4.2072 

1.9282 

115.5 

0.078125 

5.98 

-4.1966 

1.9215 

249.4 

0.093750 

4.98 

-4.1786 

1.9138 

368.6 


TABLE II. Energy (Eq/N), free energy correction 
(Esw/A’/cbT = ®P^^ tran¬ 

sition temperature (T^rlns) Q-s defined in Eq. (18) 
for different wave numbers (q) and corresponding 
wavelengths (A) for the FM-SS I transition, cal¬ 
culated for a lattice of iV = 128 x 64 atoms with 
periodic boundary conditions. 


q (tI^) ^ (nm) ^trans (K) T/SL (K) To (K) 


0.000000 

00 

0.0 

0.0 

201.8 

0.031250 

14.94 

49.4 

38.4 

223.9 

0.062500 

7.47 

88.5 

91.1 

261.7 


TABLE III. Transition temperatures for different 
wave numbers (q) and corresponding wavelengths 
(A), calculated for a lattice size of iV = 64x32 with 
periodic boundary conditions. Ttrans and TSL in- 
dicate the temperature where the SS I spiral with 
the given wave number becomes the global mini¬ 
mum of the free energy derived from the perturba¬ 
tion theory, Eq. (19), and for free spin waves, Eq. 

respectively. Tc is the temperature where the 
state becomes unstable according to perturbation 
theory. 


ity and becomes paramagnetic. This gives 
an upper bound for the transition tempera¬ 
tures, Ttrans- Ttic rcsults are summarized in 
Table [TTI| It is worth noting that although the 
ferromagnetic state remains metastable for a 
wide temperature range in the SS I phase, 
there is a temperature region where only the 
spiral state is stable and the ferromagnetic 
state becomes paramagnetic, in agreement 
with the prediction of Ref. ES] Including the 
perturbative correction also modifies the tran¬ 
sition temperature Ttrans compared to the non¬ 
interacting case. The transition temperature 
from the ferromagnetic state to the first spiral 
state is significantly increased for the interact¬ 
ing case, as can be inferred from Eig. and 
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FIG. 8. (color online) Temperature dependence 
of free energy differences AF/N between differ¬ 
ent states, including the ferromagnetic state and 
spin spirals at different wave numbers. The line 
at AF/N = 0 is a guide to the eye, identifying the 
transition temperatures. The differences obtained 
with perturbation theory, Eq. (19), are compared 
to the linear functions of the free spin wave theory, 
Eq. (17), for a lattice size of N = 64 x 32. The 
wave numbers are given in units of 


Table iml Interestingly, the transition temper¬ 
ature between the spin spiral states with dif¬ 
ferent wave vectors is hardly affected by the 
perturbation correction. Note that the 
transition temperatures are slightly different 
in Tables El and uni because of the different 
lattice sizes used in the calculations. The rea¬ 
son for this is that the lattice size influences 
not only the allowed q values, but also the spin 
wave energies. 


Fig. shows the results of Monte Carlo sim¬ 
ulations for the same model system. As clear 
from Fig. [^a), the simulation results are in 
good agreement with the spin wave calcula¬ 
tions: starting from a ferromagnetic ground 
state, the system will turn into a spiral state 
with gradually increasing wave vector until 
the temperature becomes high enough to re¬ 
move all kinds of magnetic order from the 
system. The values for /ceTtrans are some¬ 
what inaccurate (compare Tables nil with 
Fig.§, mainly because the transitions appar¬ 
ently show hysteresis. The lower wave vector 
states will remain metastable at higher tem¬ 
peratures than the point where the free en¬ 
ergy minimum moves to a different wave vec¬ 



T(K) 

FIG. 9. (color online) The order parameter de¬ 
fined in Eq. 0 for different wave numbers ob¬ 
tained from Monte Garlo simulations as a function 
of (a) increasing and (b) decreasing temperature, 
describing the EM-SS I transition, for a lattice size 
of A = 128 X 64. The wave numbers are given in 
units of ^. 

V 2a 


tor (see Fig.|^. This is even more pronounced 
in Fig. ib), where the simulation was per¬ 
formed for decreasing temperature, starting 
from a random initial state. Although the 
q = 0.062500 state is not the ground state, 
the system freezes into this metastable state 
in this case. On the other hand, the tran¬ 
sition point to the paramagnetic state is 
well approximated by the perturbation theory: 
for the q = 0.062500 spiral state, it predicts 
Tc = 261.7 K, while the critical temperature 
from the simulation is around 220 K. For com¬ 
parison, the random phase approximatiorP^ 
gives Tc = 271.1 K for the critical tempera¬ 
ture of the ferromagnetic state. The same kind 
of transition was obtained using the ab initio 
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^ ^ ^ (mRyd) (K) 

0.593750 0.79 -2.9894 1.4779 0.0 

0.156250 2.99 -2.9736 1.4168 40.8 


TABLE IV. Energy (Eq/V), free energy correc¬ 
tion {Fsw/NkBT = Incjfc/V) per spin and 
transition temperature (T/^^lns) values as in Ta¬ 
ble]^ for different wave numbers (or wavelengths) 
for the SS II-SS I transition, for a lattice size of 
V = 128 X 64 with periodic boundary conditions. 


coupling coefficients instead of the model pa¬ 
rameters, compare Fig. [^b) with Fig. [^a). 


B. The SS II-SS I transition 



T(K) 

FIG. 10. The order parameter for different wave 
numbers as a function of temperature obtained 
from Monte Carlo simulations for the SS II-SS I 
transition, for a lattice size of V = 128 x 64. The 
wave numbers are given in units of 


The SS II-SS I transition can be examined 
using the same methods as in the previous 
case. The main difference is that the energy 
of the spin spiral must have two different min¬ 
ima, both corresponding to spiral orderings, 
that is qi,q 2 ^ 0 (see Fig. [^. This requires 
at least four different coupling coefficients in 
the spin model along the x axis, illus¬ 
trated in Fig. For the model calculations 


we chose Ji = —2.0mRyd, J 3 = 2.58 mRyd, 
J 7 = —1.0 mRyd and Jn = 0.8 mRyd, which 
could reproduce the shape of the curves in 
Fig. with a slightly lower minimum at high 
wave number q = 0.593750 and a somewhat 
higher one at g = 0.156250. We also consid¬ 
ered the same Dzyaloshinsky-Moriya interac¬ 
tion between the nearest neighbors as in the 
previous case, Di = 0.4 mRyd, since this is 
necessary to stabilize both spiral states at zero 
temperature, see Eq. (14) and the subsequent 
discussion. We omitted the anisotropy terms 
needed to make the ferromagnetic state ener¬ 
getically favorable in Sec. img since they are 
irrelevant for the current discussion. The en¬ 
ergies and free energy corrections are given in 
Table [TV| for lattice size of = 128 x 64. 

The spin wave calculations indicate that 
starting from a high wave vector ground state, 
the system may indeed switch to a low wave 
vector ordering. This is in agreement with the 
Monte Carlo simulation results with the same 
parameter set, shown in Fig. [M as well as 


simulations performed with the ab initio cou¬ 
pling coefficients, see Fig. |^c). The spin wave 
expansion again underestimates the transition 
temperature as in the case of the FM-SS I tran¬ 
sition. 


C. The SS III-SS II transition 

The third type of transition found in the 
Fe monolayer on Ta(llO) surface corresponds 
to the case when the wave vector of the spi¬ 
ral remains fixed, but the spiral normal vector 
rotates from the y axis (the cycloidal state) 
towards a direction in the xy plane. For mod¬ 
elling this transition we supposed that the 
wave vector q^ of the spiral state is determined 
by the isotropic exchange couplings, while 
the Dzyaloshinsky-Moriya interaction and the 
anisotropy terms were taken into account as 
a perturbation. For the anisotropy we chose 
Kx > 0 and Kz = 0 , since ab initio calcula¬ 
tions indicated that at 15% relaxation the fer¬ 
romagnetic states along the y and 2 ; axes have 
almost the same energy, while the x axis is a 
hard axis (see Fig. ^b)). The angle between 
the xz plane and the plane of the spin spiral 
will be denoted by ip. In this case, the energy 
contribution per spin from the Dzyaloshinsky- 
Moriya interaction and the anisotropy terms 
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can be expressed as 

^ = - 5-0 (go) cos if + cos^ <^.(20) 

Differentiating ( [2Q| ) with respect to ip leads 
to the stationary points 


simp^^^ = 0 , ( 21 ) 


Substituting the solutions into (20) gives 


N 

A£’(2) 

N 


—^ 2 ^^ (9o) + 2 ’^® 


8K, 


(23) 

(24) 


implying that whenever the second stationary 
point exists, 


i-0 (go) 


2K, 


< 1 , 


(25) 


it will correspond to the energy minimum. 
This describes the rotation of the spiral 
normal vector away from the y axis when 
the Dzyaloshinsky-Moriya interaction is weak 
compared to the anisotropy. Calculating the 
spin wave spectrum reveals that only one of 
the states is stable for any value of D and 
therefore the spin wave expansion is not suit¬ 
able for describing this type of transition. 

For the present simulations the exchange pa¬ 
rameters Ji = — 2.0mRyd, J 3 = S.OmRydand 
Jj = — l.OmRyd were chosen which lead to a 
spin spiral along the x axis with a wave num¬ 
ber q = 0.546875 (A = 0.85 nm). We took 
Di = 0.05 mRyd between the nearest neigh¬ 
bors and Kx = 0.2 mRyd, and found that 
these values did not influence the shape of the 
spiral considerably, but confined the spins to 
a plane with a normal vector lying in the xy 
plane, as shown in Fig. Sd)- We also used 
a ferromagnetic coupling between the neigh¬ 
bors in the y direction, J 2 = —2.0 mRyd, 
which does not influence the spiral state but 
removes the possible domain walls from the 
system along the v ax is. These domain walls 
occur because Eq. (22) has two solutions 


with the same energy, therefore if the spins are 



T(K) 

FIG. 11. (color online) Different components of 
the order parameter Eq. ([^ as a function of tem¬ 
perature, for the SS III-SS II transition. The wave 
number of the spin spiral was q = 0.546875^^ 
and a lattice size of A/" = 128 x 64 was used. 


weakly coupled along the y direction, and 
—ip^‘^^ domains may be simultaneously present 
in the system. 

The SS III-SS II transition is shown in 
Fig. [TTJ in agreement with the simulations 
performed with ab initio coupling coefficients, 
see Eig. id). By increasing the temperature, 
the plane of the normal vector of the spiral 
rotates towards the y axis, which is the one 
preferred by the Dzyaloshinsky-Moriya inter¬ 
action over the x direction preferred by the 
anisotropy. This indicates that with increas¬ 
ing temperature the magnitude of the effective 
Dzyaloshinsky-Moriya contribution to the free 
energy decreases slower than the anisotropy 
contribution. 


IV. SUMMARY AND CONCLUSIONS 

We examined the phase diagram of an Ee 
monolayer on the (110) surfaces of W and Ta 
as a function of the relaxation of the Ee layer 
and the temperature. We used the relativistic 
screened Korringa-Kohn-Rostoker method to 
determine the single-particle potential of the 
systems within the local density approxima¬ 
tion of density functional theory. In terms of 
the relativistic torque method, we calculated 
the tensorial coupling coefficients which ap¬ 
pear in the generalized Heisenberg model de- 
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scribing the spin system, Eq. 0 . Based on 
this spin model, we determined the magnetic 
ground state from spin dynamics simulations, 
and we performed Monte Carlo simulations to 
explore the magnetic phase transitions at fi¬ 
nite temperature. 

In case of W substrate the obtained mag¬ 
netic moments and ground states were in good 
agreement with previous calculation^®^ and 
with experiment s.l^ The ground state was 
ferromagnetic with an easy axis along the 
[110] direction for relaxations smaller than 
15%, including the experimentally and theo¬ 
retically determined relaxation values around 
12% — 13%. For larger relaxations, the sys¬ 
tem remained ferromagnetic, but the easy axis 
turned into the out-of-plane [110] direction. 
For fixed relaxations, we found no thermally 
induced transition between these two states. 

In case of Ta substrate four different phases 
were identified in the considered relaxation 
range, see Fig.l^b), with transitions occurring 
at 10.5%, 13.8% and 14.5% relaxation values. 
At low relaxations the ground state was ferro¬ 
magnetic with an easy axis along [110]. The 
next two phases, denoted by SS I and SS II, 
correspond to cycloidal spin spirals with wave 
vectors along the [110] direction and normal 
vector along the [001] axis, the SS II state hav¬ 
ing a significantly larger wave number. The 
SS II and SS III spin spirals had similar wave 
vectors but the normal vector of the spiral left 
the [001] axis in the SS III state. 

Choosing the relaxation value close to one 
of the transition points, different types of 
transitions were obtained between these states 
at finite temperature. These possible phase 
transitions were described theoretically using 
spin wave expansion and compared to Monte 
Carlo simulations performed on model sys¬ 
tems. We found that starting from a ferro¬ 
magnetic ground state, the system may turn 
into a spin spiral state at finite temperature 
before becoming paramagnetic. Although the 
appearance of the spin spiral state as a conse¬ 
quence of the Dzyaloshinsky-Moriya interac¬ 
tion is a well-known effect in two-dimensional 
systems such as Mn monolayeJ^ or Fe double- 
layei^ on W(llO), there was no such tran¬ 
sition observed as a function of temperature. 
However, ab initio calculation^!®^ indicated 


a ferromagnetic ground state for Fe double¬ 
layer on W(llO), suggesting that this system 
is probably very close to such a ferromagnetic- 
spin spiral transition. We have also shown 
that the high wave vector SS II state may turn 
into the low wave vector SS I spiral by increas¬ 
ing the temperature, while in the case of the 
SS III-SS II transition the normal vector of 
the spin spiral rotated from a general in-plane 
direction towards the [001] direction. 

For all three phase transitions, the simula¬ 
tions performed on model systems and using 
the ab initio coupling coefficients gave results 
which were in agreement with the predictions 
based on spin wave expansion. Compared to 
the Monte Carlo simulations, the spin wave ex¬ 
pansion gave good approximations for the tem¬ 
perature Tc where any magnetic order disap¬ 
pears and somewhat underestimated the tran¬ 
sition temperature Ttrans between the ordered 
states. The latter difference is also a conse¬ 
quence of the metastability of the states, indi¬ 
cating that conventional Monte Carlo simula¬ 
tions are not well-suited for finding the actual 
transition temperature. 

Given the wide variety of possible ground 
states in a relatively narrow range of relax¬ 
ations, our present work might motivate ex¬ 
periments to determine the actual magnetic 
ground state of an Fe monolayer on Ta(llO). 
It may even be possible to find one of the 
thermally induced transitions described here. 
On the other hand, the spin wave expansion 
method may also be applied for the finite- 
temperature description of other stable equi¬ 
librium configurations such as the skyrmion 
lattice structure found in ultrathin magnetic 
films PH 
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